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Abstract 

These notes are an introduction to a few selected theoretical ideas in the field of quantum 
spin liquids: classical zero modes and breakdown of the 1/S expansion, the Lieb-Schultz- 
Mattis-Hastings theorem and Oshikawa's argument, the short-ranged resonating valence-bond 
picture, large- A*' limit (Schwinger bosons) and Z2 gauge theory. 
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1 Introduction: band and Mott insulators 

Depending on the context (experiments, theory, simulations,...), "Quantum spin liquid" is 
sometimes used with rather different meanings. But let us start with a first simple definition: 
the ground state of a lattice quantum spin model is said to be a quantum spin liquid (QSL) 
if it spontaneously breaks no symmetry. According to this first definition, a QSL is realized if 
the spins fail to develop any kind of long range order at zero temperature (T = 0) (hence the 
word "liquid", as opposed to solids which are ordered and break some symmetries). Of course, 
this first definition raises a number of questions: Does this define new distinct states of matter 

* Lectures given at the Les Houches summer school on "Exact Methods in Low-dimensional Statistical Physics 
and Quantum Computing" (July 2008). 
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? Do QSL have some interesting properties ? Are there some experimental examples ? To 
answer these questions, it is useful to go back to the origin of magnetism in insulators. 

Generally speaking, there are two kinds of insulators: band insulators, and Mott insulators. 
The first ones can be qualitatively understood from the limit of non-interacting (or weakly 
interacting) electrons. Consider for instance a periodic lattice]^ with an even number n of 
sites per unit cell, with an average electron density of one electron per site (so-called half 
filling). The Hamiltonian describing how the electrons hop from sites to sites looks like 

Hk = ~iX](ij)CT=Ti {^iir'^^'^ ~^ ^■'^ ' "'^here only first neighbor hopping is considered for 
simplicity. H can be diagonalized in Fourier space and gives n dispersing bands. The ground 
state is just the Fermi sea obtained by filling the lowest energy states. Since the density 
is one electron per site, the n/2 lowest energy bands are completely filled (one up and one 
down electron for in each single particle state). Assuming that the band n/2 -f 1 is separated 
by a gap A in energy from the n/2 lower bands, all the excitations are gapped and, at 
temperatures smaller than the gap, there is no charge carrier to carry an electric current. 
This is the well known picture for a band insulator: there are no low energy charge degrees 
freedom, no magnetic (spin) degrees of freedom, the ground state (Fermi sea) is unique and 
breaks no symmetry. To get some interesting QSL, we should instead look at Mott insulators. 
There, the number of sites per unit cell is odd and the non-interacting limit is unable to 
give the correct insulating behavior (at least one band is partially filled, hence with low 
energy charge excitations). It is more useful to look at the system in the opposite limit 
of very large electron-electron repulsion, as with the large U limit of the Hubbard model: 
H = Hk + U '^li '^iT ^li'^il ■ At U = oo and t = (still at half filling), the ground state 
is highly degenerate (= 2^, where V is the total number of sites) since any state with one 
electron per site is a ground state, whatever the spins orientations. To describe how this 
degeneracy is lifted at weak but finite t/U, a second order perturbation has to be computed.^ 
The result is an effective Hamiltonian acting in the subspace spin configurations, and takes 
the form of a quantum spin-^ Heisenberg model: 

H = - ^ JrjSr ■ Sj (1) 



where Jij = t^j/U involves the hopping amplitude tij between sites i and j and measures the 
strength of the antiferromagnetic ( AF) interaction between the (electron) spins Si and Sj |^ 

Although the model of Eq. |^ is in general a complicated quantum many body problem 
with very few exact results,^ its ground state and low energy properties are qualitatively well 
understood in many cases. In particular, the ground state can be antiferromagnetically ordered 
(also called Neel state). Such state can be approached from a semi classical point of view 
described in Sec. ^ the spins point well defined directions and form a regular structure. Most 
of the Mott insulators studied experimentally belong to this family. The simplest example 
is the nearest neighbor Heisenberg model on bipartite lattices such as the square, cubic or 
hexagonal lattices. There, on average, all the spins of the sublattice A point in direction 
-I- So (spontaneous symmetry breaking of the SU{2) rotation symmetry), and all the spins of 
sublattice B point in direction —So- The difference with a classical spin configuration is that 
the magnetization of one sublattice (it is the order parameter for a Neel state) is reduced by 
the quantum zero-point fiuctuations of the spins, even at T = 0. Such ordered states are not 
QSL (they might instead be called spin "solids") since they break the rotation symmetry. 

The main question addressed in these notes is the fate of the ground state of Eq. |l|whe n the 
lattice and the interactions Jij are such that the spins fail to develop any such Neel ordered 
state. A state without any order is not necessarily interesting from a theoretical point of 



^We use a tight binding model where the sohd is modeled by one state per site, neglecting (or, more precisely, 
integrated out) filled orbitals or high energy empty states. 

^At first order in t, a single electron hopping inevitably leads to a doubly occupied site. 

^In real materials, there are often tens or hundreds of electron per unit cell, several ions and many atomic 
orbitals. Although the description of the magnetic properties in terms of lattice spin models if often very accurate, 
the spin-spin interactions is often more complicated than this antiferromagnetic Heisenberg model. It is quite 
frequent that some interactions violate the SU{2) symmetry of the Heisenberg model, due to spin-orbit couplings 
in a crystalline environment. In these notes, we focus on models with an SU{2) symmetry. 

■^In these notes, we focus here on dimension D > 1, but much more is known about one-dimensional (ID) spin 
chains. 
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view. For instance, a spin system at very high temperature is completely disordered and does 
not have any rich structure. As we will see, the situation in Mott insulators at T = is 
completely different. A first hint that Mott QSL host some interesting topological properties 
will be discussed in Sec. ^ (Lieb-Shultz-Mattis ||l| Hastings |^ theorem). A concrete (but 
qualitative) picture for QSL wave functions is given in Sec. H, in terms of short range valence 
bond configurations and deconfined spinous (magnetic excitations carrying a spin ^). Finally, 
Sec. ^, presents a formalism which puts some of the ideas above on firmer grounds. It is 
based on a large-A'' generalization of the Heisenberg models (5(7(2) Sp{N)) which allow 
to describe some gapped QSL and to establish a connection by topologically ordered state of 
matter, such as the ground state of Kiatev's toric code j^. 



2 Some materials without magnetic order at T = 

There are many magnetic insulators that do order at T = 0.^ For instance, the magnetic, 
properties of many compounds are described by ID spin chains of spin ladder Hamiltonians. 
Thanks to the Mermin- Wagner theorem and the reduced dimensionality, these system cannot 
develop long range spin-spin correlations, even at T = oj^ They certainly deserves to be called 
QSL and represent a very rich field of activity. In these notes we will instead focus on QSL 
in D > 1 systems, where our present understanding is less complete. 

CaV4 09 is the first Heisenberg system in D > 1 where the magnetic excitations were 
experimentally shown to be gapped, in 1995 |^. This compound can be modeled by an an- 
tiferromagnetic spin-i Heisenberg model on a depleted square lattice where one site out of 
five is missing (Fig. |]|). The remaining sites correspond to the locations of the Vanadium 
ions, which carry the magnetically active electrons (spins). The magnetic interactions Jij 
turned out to be significant not only between nearest neighbors, but also between second 
nearest neighbors (the electron hops through oxygen orbitals, which have a complex geome- 
try). Through magnetic susceptibility measurements, it was shown that the ground state is 
a rotationally invariant spin singlet, thus excluding any Neel ordering. This QSL behavior 
can be understood by taking a limit where only the strongest Jij are kept, and the other are 
set to zero. It turns out these strongest couplings are between second-nearest neighbors, and 
form a set of decoupled four-site plaquettes (of area \/2 x \/2 and surrounding a missing site). 
Since the ground state of such a four-site Heisenberg cluster is a unique singlet S = state, 
separated by a gap from other states, the model is trivially a gapped and without any broken 
symmetry in this limit. But this is not the kind of QSL we want to focus on here, since it 
can be adiabatically transformed into a band insulator. Switching off the electron-electron 
interactions would make the system metallic, but one can proceed in a different way. Starting 
with realistic values of the Jij , the inter plaquette couplings are gradually turned off. Doing 
so, one can check (numerically for instance) that the spin gap does not close and no (quantum) 
phase transition in encountered. Then, in this systems of decoupled four-electron cluster, the 
Hubbard repulsion U can be switched to zero, without causing any phase transition. The final 
model is evidently a band insulator and smoothly connected to the initial Heisenberg model. 

Since then, numerous 2D and 3D (Heisenberg) magnetic systems with an even number 
of spin-i per unit cell have been found to be gapped. To our knowledge, their ground state 
can be qualitatively understood from a limit of weakly coupled clusters in all cases and can 
therefore be "classified" as band insulators (as CaV4 09 above). Some of them can be very 
interesting for different reasons,^ but their ground states are not fundamentally new states of 
matter. 

In the recent years, experimentalists have also uncovered a number of materials which are 
well described by 2D Heisenberg models with an odd number of spin-| per crystal unit cell,^ 
and which do not develop any Neel order when T ^ 0. Some examples are the Herbertsmithite 

'Some order at a temperature with is very small compared the typical energy scale of the Heisenberg spin-spin 
interactions. This is often due to perturbations that are not included in the simplest Heisenberg model description. 

^Due to some residual 3D couplings, there can be a finite temperature phase transition to an ordered state at 
very low temperature. 

'^For instance: TlCuCla jij is coupled dimer system with a Bose-Einstein condensation of magnetic excitations 
in presence of an external magnetic field, and SrCus (603)2 |p|| has a magnetization curve with quantized plateaus. 
*In such case, the absence of long range order cannot be attributed to some band insulator physics. 
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Figure 1: Depleted square lattice model for the magnetic properties of CaV40g. The different 
exchange energies are shown by different types of line. The strongest J correspond to the fat Hues 
forming the large tilted square plaquettes. 



(ZnCu3(OH)6Cl2)| §1 and Volborthite (Cu3V207(OH)2 •2H2 0)||lminerals (both with a 
kagome lattice geometry), triangular based organics materials |9|, [L0|, or triangular atomic 
layers of He^ adsorbed onto graphite |0 (there the spin is not electronic, but nuclear). It turns 
out that all these systems seem to have gapless magnetic excitations and a complete theoretical 
understanding of these system is still lacking. The present theories for gapless QSL are rather 
elaborate [ p^ and many questions remain open (stability, nature of the excitations, correlation 
exponents, etc.). However, as we will see, gapped QSL are simplest from a theoretical point 
of view. Intriguingly, to our knowledge, no gapped QSL has been discovered so far in nature, 
although many spin models do have gapped QSL ground states. 

3 Spin wave theory, zero modes and breakdown of 
the 1/S expansion 

To understand why an AF Heisenberg spin model can fail to order at zero temperature, is is 
useful to briefly review the standard approach to Neel phases: the semi-classical 1/S' spin- wave 
expansion jl3||. This approach i) starts from a classical spin configuration which minimizes the 
Heisenberg interaction ii) assumes that the quantum deviations from this ordered direction 
are small iii) treats this deviations as collection of harmonic oscillators (the leading term in a 
1/S expansion). In this approximation the Hamiltonian is written using boson creation and 
annihilation operators, is quadratic, and can be diagonalized by a Bogoliubov transformation. 
One can then check a posteriori if the spin deviations are indeed small. If it is not the case, 
we have a strong indication that the magnetic long range order is in fact "destroyed" by the 
quantum fluctuations, thus opening a route for a QSL ground state. 

3.1 Holstein-PrimakofF representation 

The starting point is the representation of the spin operators using Holstein-Primakoff jlj] 
bosons 

= S- a\a, , S+ = sj2S^alar Ui , S' ^ a\^2S~ala„ (2) 

from which on can check that the commutation relations [Sf, Sf ] = ie'^'^^ Sf and Si = 5(5* -fl) 
are satisfied (using [aijaj] ~ 1). 

Let {zi} be a classical ground state of Eq. |l|, minimizing E = ^ YLij JijSi ■ zj with = 1. 
These directions can be used as local quantization axes: we use Eq. in a local (orthogonal) 
frame {xi,yi,Zi = Xi A i/i) adapted to the classical ground state. Under the assumption that 
Si shows small deviations from the classical vector Szi, the typical number (a^a) of Holstein- 
Primakoff bosons should be small compared to S. We can therefore simplify S~^ (and S~ ) in 



^Although the spin-spin interaction strength is of the order of J ~ 200 K, no order has be found down to 50mK. 



Eq. g by keeping only \/2S' in the square roots, to obtain ||l 
where 



-2 t 1 



and Zi ■ Tii — 0. 
Replacing Eq. |3| in the Hamiltonian gives 

ij ij 



(3) 



(4) 

(5) 
(6) 



(7) 



The first term is a constant, proportional to the classical energy Eo- The two other terms, 
proportional to S, are quadratic in the boson operators and describe the spin fluctuations 
as a set of coupled harmonic oscillators.^ The positions qi = -^{ui + a\) and momenta 

Pi = ~ 0,1) operators of these oscillators can be conveniently grouped into a column 

vector of size 2N (N is the total number of spins): 



so that H becomes 



H - 



Pi 

PN 



(S+ifi5o + |v*MV, 



where A4 is a 2N x 2N matrix given by 



M 



J' 



jyy _ jz 



and the N X N matrices .r^^ , -F^ , .F^ and are deflned by: 



'-'ij — 'Jij ^i ' X 



and Jif — Sij ^ JikZi ■ 



Zk- 



(8) 



(9) 



(10) 



(11) 
(12) 



3.2 Bogoliubov transformation 

Diagonalizing H amounts to find bosonic creation operators and corresponding energies 
LOa > such that H = cOa {baba + I) (up to a constant). A necessary condition is that the 
operator 6j, and ba are "eigenoperators" of the commutation with H, for the eigenvalues ujc, 
and —uja respectively: [_ff, bj,] — ujabl, and [H, We thus seek the eigenvectors 

of the action of [H, •] in the space of linear combinations of qi and pj . The commutators of 
H (Eq. m with the operators q and p are simple to obtain using [qi,qj] = ~ and 



^"Due to the fact that {ii} minimizes the classical energy, -hjSj is perpendicular to Zi and thus orthogonal 
to TTi, and there is no term linear in tt. 
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[litPj] = "iSij. For an arbitrary linear combinations of the Qi qnd Pi parametrized by the 
complex numbers xi, - ■ ■ , X2n the result is 



[H,Xiqi + xnQn + XN+iVi H X2nPn] 

= yigi + yjvgiv + j/iv+ipi H VznPn 



with the coefficients yi , ■ ■ • , j/2jv given by 





^iS M 


" 1 " 




Xl 


. l/2]V . 




_ -1 




X2N 



(13) 



(14) 



where 1 is the N x N identity matrix. So, finding the operators (spin-wave creation 
operators) amounts to find the eigenvectors of the "commutation matrix" C = iM ^ 

But C is not symmetric and cannot always be fully diagonalized (contrary to TVf). It can 
be shown that if all the eigenvalues of M were strictly positive, C could be diagonalized, its 
eigenvalues would be real and come in pairs —a), a) J"] 

However, Ai does have some zero eigenvalues. The matrix Ai is not specific to the 
quantum spin problem. The quadratic form describing the classical energy variation for a small 
perturbation around the chosen classical ground state {zi}, is described by the same matrix 
In particular, if the classical ground state admits some zero energy (infinitesimal) spin 
rotations, M posses some eigenvector for the eigenvalue 0. Because global rotations should 
not change the energy, A4 has at least two zero eigenvalues. Still, these global rotations do not 
cause difficulties in diagonalizing the spin-wave Hamiltonian, they just correspond to some 
uJct — (the associated collective coordinate Q and conjugate momentum P simply do not 
appear in H). 



3.3 Zero modes on the kagome lattice 

However, some Heisenberg models admit classical zero modes (hence zero eigenvalues in A^) 
which do not correspond to global rotations. As an example, consider the Heisenberg model 
on the kagome lattice ||li[ (for another classic example, the J1-J2 model on the square lattice, 
see llill ). Any classical spin configuration such that the sum Zi + Zj + Zk vanishes on each 
triangle {ijk) minimizes the classical energy. Among the numerous ways to achieve these 
conditions, are the planar ground states, where all the spins lie in the same plane. In such 
a state, the spins take only three possible directions, a, b and c at 120 degrees from each 
other. On the kagome lattice, there is an exponential number of ways to assign these three 
orientations such that the same letter is never found twice on the same triangle (three-coloring 
problem, see Fig. Now, choose one of these "a&c" states, and find a closed loop of the type 
ababab • • • . Because of the three-coloring rule, the spins which are neighbors of this loop all 
point in the c direction. Now, we can rotate the spins of the loop about the caxis by any angle. 
This transforms the planar ground state into another (non planar) ground state, without any 
energy cost. So, for a generic planar ground state, we get as many zero modes (in Ai) as 
closed loops with two alternating "colors". This number typically grows like the number of 
sites in the system. 

What are the consequences of such classical zero modes for the quantum problem ? As 
explained previously, the operators describing the two transverse directions along which the 
spins can deviate from the Zi axis obey the same commutation rules (at leading order in 
the l/S expansion) as the position q and momentum p of an harmonic oscillator. In the 
case of the kagome loop modes discussed above, the energy is zero in one direction (rotation 



^^Let P be an orthogonal matrix which diagonalizes symmetric M: M = P ^AP, where A is a diagonal matrix 
and PP* = 1. If the eigenvalues of M are strictly positive, K = P^^\/\ is invertible and M = KK^. We 

" Then, C = K~^CK = iS K* aK is Hermitian (since a is real 



write C = iS KK a, where u 



-1 



antisymmetric, and K is real). C can therefore be diagonalized and its spectrum is real. Since C and C have the 
same spectrum, C can also be diagonalized and has real eigenvalues. Finally, we use C* = —C. Since C and C* 
should have the sanie spectrum, the eigenvalues of C go in pairs — aj,aj. 

12 The Eqs. U and M also hold if TTi is a classical spin deviation of length vrj <C 1. 
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Figure 2: Classical planar ground state on the kagome lattice. The loops where the spins alternate 
between the a and b directions are marked with dashed lines, they host independent zero modes 
(by rotation around the c axis) . 



about the c direction), and quadratic in the other direction. Using the associated collective 
coordinate P and Q, we expect the Hamiltonian to be proportional to H = ^{P^ +lj^Q'^) 
with u! = 0, since there is no classical energy cost for spin deviation in the direction Q. The 

"0 u" ~ 



corresponding commutation matrix is C = i 



-1 



and cannot be diagonalized when 



u! = 0, as anticipated. In general, each such local zero mode will lead to an irreducible 2x2 
Jordan block of this kind.[j The ground state |0) of the oscillator is simple to obtain and 
corresponds to a zero point motion of the coordinate Q which diverges when uj ^ (no 
restoring force, like for a free particle) : {0|(5^|0) = 

As long as the the number of such zero modes is finite in the thermodynamic limit (this 
is the case when the classical ground state has no special degeneracy, beyond those implied 
by global rotations), the divergences above have a zero measure and do not cause divergences 
in the number of bosons (0|ajai|0),[^ which measures the strength of the deviations from 
the classical state. In such a case, the Neel ordered state is stable with respect to quantum 
fluctuations, at least for large enough On the other hand, if the number of such modes 
grow like A*', the average number of bosons diverge and the spin-wave expansion breaks down 
(the initial assumption that {0|ajai|0) is flnite and small compared to S cannot be satisfled). 

At this point, a route to obtain a QSL appears to look for a lattice where the classical 
model has a suflEicient number of "soft" modes, so that the zero point motion of the spins restore 
the rotation invariance and destroy the long range spin spin correlations. This condition is 
realized on the kagome lattice, where indeed all numerical studies concluded to the absence 
of Neel order in this system (at least for S — |). However, the semi classical spin wave 
theory described here breaks down. As discussed in the next sections, QSL states in Mott 
insulators possess some internal topological properties which are missed by the simple picture 
of a "disordered" state which would just be the quantum analog of a high temperature phase. 



^^The general theory for possible Jordan forms of C (size and nature of the irreducible blocks) is in fact a result 
of classical mechanics, found by Williamson and exposed in m^- 

^* {0| aja^ 1 0) = ^(0|p| +qj ~ 1|0) can be computed by expressing qt and Pi in terms of bJ, and ba, or in terms of the 
new position and momenta Qa = '^{^a + ba) and Pa = "^^(''q — 6a) Concentrating on the term (0|q^|0), qi is a 

linear combination of the type qi = J2a=i ^aQa~^J2^=i '"^^Pp^ and v are related to the eigenvectors of C). From 
the fact that |0) is the vacuum of the ba bosons, we have (0|PiPj|0) = {0\QiQj\0) if i ^ j, and {0\PiQj+QjPi\0) = 
Vi, j. Then the square of the spin deviation at site i (here the Xi component) is a linear combination of the zero 
point fluctuations of the normal harmonic oscillators (0|i}||0) = l^aC^D^CIQalO) + Z]a ('^D^ |0) • Assuming 
a regular behavior of the coefficient (m^)^ and (Va)^, (0|g^|0) is typically the sum of terms proportional to ~ l/wa 
when the mode frequency uia is small. 

^^This does not imply that the order should persists down to S = ^. 
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Figure 3: A lattice model which is translation invariant and periodic in the x direction can be 
viewed as a ring. The interactions Jy , indicated by dashed lines, are invariant in the x direction 
but otherwise arbitrary. In this example, each cross section has C — 3 sites. 



4 Lieb-Schultz-Mattis theorem, and Hastings's exten- 
sion to D > 1: ground state degeneracy in gapped spin 
liquids 

The Lieb-Schultz-Mattis theorem |Q] was originally derived for spin chains and spin ladders 

f |l9| and was recently extended to higher dimensions in an important work by Hastings 
(see also jioj for an intuitive topological argument valid in any dimension, and jilj for 
a mathematically rigorous proof). It applies to spin Hamiltonians which are translation 
invariant in one direction (say x), have a conserved magnetization Sfot = , and short 

range interactions. In addition, the model should have periodic boundary conditions in the x 
direction. Although more general interactions can easily be considered,rj we concentrate for 
simplicity on spin-S Heisenberg models, written as in Eq. |l] (with Jij — Ji+x,j+x to respect 
the translation invariance). 

Following (iot , we define the cross section as all the sites sitting at a given value of x. By 
translation invariance, all cross sections are equivalent and contain C sites (Fig. H). In a spin 
chain, each cross section contains a single site. In an n— leg spin ladder, C — n sites. In a 
square lattice, C = Ly. On a D-dimensional lattice with n sites per unit cell, C = nL^~^ , 
etc. We note Lx the system length in the x direction, and therefore CLx is the total number 
of sites. Finally we define = {Q\Sia\\Q) as the ground state magnetization per site. 

The theorem says that if (7(5 + m^) is not an integer, the ground state is either degenerate, 
or the spectrum has gapless excitations in the thermodynamic limit. In other words, if 
C{S -\- m") ^ Z the system cannot have a unique ground state and a finite gap to excited 
states in the thermodynamic limit. Although the proof in ID |l|| and Oshikawa's topological 
argument ||2^] (Sec. 0) are relatively simple, the proof appears quite involved for D > 1, and 
will not be discussed here. 

What is the relation between the LSMH theorem and QSL ? In most AF Heisenberg 
models on a finite-size lattice, |0) is a singlet and — 0. If we focus on the case S = ^, 
the theorem forbids a single ground state and a gap when C is odd. In particular, if the 
lattice is two dimensional and describes a Mott insulator, the unit cell has an odd number 
n of sites and any odd Ly can be chosen to get and odd C = nLy (note that the total 
number of sites is still even if Lx is even). If we assume that a gapped QSL is realized 
(for an example which fits in the LSMH conditions, see for instance |22|]), its ground state 
must be degenerate (with periodic boundary conditions). Usually, ground state degeneracies 
are the signature of some spontaneous symmetry breaking. However, by definition, a QSL 
respect all lattice symmetries. The degeneracy imposed by the LSMH theorem cannot be 
understood from this conventional point of view and is a hint that (gapped) QSL wave function 
possess some interesting topolo gica l properties, which correspond to the notion of "topological 
order" introduced by Wen |p3|, |24|| for spin systems and Wen and Niu | p5| | in the context of 
the fractional quantum Hall effect. As we will briefly discuss at the end, this topological 

^^In particular, the interaction can be anisotropic: SfS^ + A{SfSj + S^S^), and an external magnetic field 
parallel to the z direction can be present. 




Figure 4: Schematic spectrum of the twisted Hamiltonian (Eq. as a function the angle 6*, in 
the case where C{S + m^) is a half integer. 



degeneracy is deeply related to the exotic nature of the elementary excitations in a 

4.1 Oshikawa's topological argument 

Oshikawa's argument is somehow related to Laughlin's argument ||2^ for the quantization 
of the transverse conductivity in the quantum Hall effect. First, a "twisted" version of the 
Hamiltonian is introduced: 



(15) 



where < Xi < is the a;-coordinate of site i. It is simple to show that the spectra of Ho 
and H2n are the same, since the unitary operator 

U = l[e^p(2in^sA (16) 

maps Ho onto H2tt: 

UHoU-^ = H2^ (17) 

(the calculation simply uses e**'^' S'+g-'^'^-" = S+e''^). 

Starting with a spectrum of Ho which is gapped, we further assume that the gap of He 
remains finite when 9 goes from to 27r.0 On can follow the ground state of He, which does 
not cross any other energy level as 9 is varied. Assuming that the ground state |0) of Ho is 
unique and using the finite gap hypothesis, it must evolve to the ground state of H2n, denoted 
|27r). Through Eq. |l|, both states are related: \2tt) = U ^|0). However, the operator U does 
not always commute with the translation operator T and may change the momentum. The 
precise relation is 

TU = UTex^p ^2i7r^^ exp (2i7rC5') . (18) 

The first phase factor, also equal to 27rCm^, comes from the shift by 2-k/Lx of the local rotation 
angles after a translation. The second phase factor corrects the 2-k jump of the rotation angle 
when passing from x = Lx ~ 1 to x — 0. This relation implies that the momentum fco of |0) 
(defined by TjO) = e'^'olO)) and the momentum k2n of |27r) = [/"^|0) are related by 

fco = fca^ + 27rC(5' + m") (19) 

But He is translation invariant (commutes with T) and the momentum of each state (quantized 
for finite L^) cannot change with 9. So |0) and \2n) have the same momentum and fco = 
k2TT [2n]. Prom Eq. |l^, we get that C{S + m^) must be a integer. 



^'^QSL have "spinons" excitations which carry a spin ^ (like an electron) but no electric charge. 

^^Hastings argument does not directly use He for a finite 9 and does not rely on this assumption. This assumption 
is however reasonable by the fact that, under an appropriate choice of gauge (frame), H = Ho and He only differ 
for the terms connecting the cross section at x = L^; — 1 to the cross section at x = (boundary terms), and are 
identical in the bulk. 
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Figure 5: A short range valence bond state on the triangular lattice. The singlet pairs are marked 
with eUipses. 



5 Anderson's short range resonating valence-bond pic- 
ture 

In the l/S expansion, it is assumed that the spins experience small fluctuations about a 
well defined direction and that spin-spin correlations are long ranged. This is of course 
incompatible with having a rotationally invariant QSL state. To gain some intuition about 
what a QSL wave function may look like, it is instructive to start from a completely opposite 
limit: a spin singlet state with extremely short range correlations. A short range valence- 
bond (VB) state is such a wave function, it is the direct product of S = states \[ij]) — 
75 (I TUj) — I iitj)) on pairs of sites : 

I VB) = I [ioii]) ® I [i2«3]) ® ! [44^5]) ® • ■ ■ I [iiv-iijv]) (20) 

where each site of the lattice appears exactly once (Fig. Such a VB state is said to be 
short range if all pairs of sites coupled in a singlet are at a distance Ir^^ — ri^_^-^ \ smaller than 
or equal to some fixed length rmax(much smaller that the lattice size). The simplest case is 
Tmax = 1, where each spin forms a singlet with one of its nearest neighbors. 

In a VB state, the spin-spin correlations are short ranged: {VB|Si • S,|VB) = if \ri^ — 
Tip+i I > '"max. For a nearest neighbor Heisenberg model on a bipartite lattice, one can compare 
the (expectation value of the) energy of a nearest neighbor VB state, with that of the simple 
two-sublattice Neel state | titi ■■■)• The VB energy is evB = ""^f P^'' site and the Neel 
one is en — ~~J^, where z is the number of nearest neighbors. If the lattice is not bipartite 
but admits a three-sublattice classical ground states (with spins pointing at 120 degrees from 
each other), the energy of a classical Neel state is eN — —Jj^- From this, we observe for 
instance that the VB energy is lower than en on the kagome lattice. More generally, this 
simple variational comparison shows that a low coordination z and frustrated interactions 
(which increase the number of sublattices in the classical ground state) tend to favor VB 
states, and thus possible QSL states. 

In fact there are many (frustrated and Heisenberg-like) toy models where some/the nearest 
neighbor VB states are exact ground states. The most famous example is the Majumdar-Gosh 
model Consider the spin-| Heisenberg chain with first- (Ji) and second- (J2) neighbor 
couplings. At Ji = 2J2 > we have 

HyiG ~ 2 Si ■ Si+\ + Si ■ Si+2 (21) 
i i 

and the (two-fold degenerate) ground states are exactly given: 

\a) = •••®|[01])®|[23])®|[46])®--- (22) 
\b) = ••■®|[12])®|[34])®|[56])®--- (23) 

The proof can be done three steps. First, the Heisenberg Hamiltonian on three site Hijk = 
Si ■ Sj + Sj ■ Sk + Sk ■ Si is written as Hijk ~ ^{Si + Sj + Sk)^ — |. In this form, proportional 
to the square of the total spin, it is clear that the eigenvalues of Hijk are ^S{S -f 1) — | 
with S = I or S = I (the only possible values of S for three spin-i). So, if the sites ijk 
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are in a. S — ^ state, they minimize exactly Hijk- Second, one expresses the Majumdar-Gosh 
Hamiltonian as 

Hmg = '^^Hi-i.i^i+i. (24) 

i 

Finally, one remarks that the dimerized states \a) and |fe) always have one singlet among the 
sites i — + 1, which are therefore in a S = 1/2 state. We conclude that \a) and \b) 
minimize all the terms in Eq. |^ and are thus ground states of Hmg- 

The Majumdar-Gosh model is the simplest model of a family o f sp in models where exact 
VB ground states can be found.|^ For instance, the Husimi cactus [ po| is a lattice constructed 
as a tree (no loops) of corner sharing triangles. Its geometry is locally similar to the kagome 
lattice but it has no closed loop (except of course for the triangles themselves). The argument 
above (writing the Hamiltonian as a sum of Hijk) directly generalizes to this case and shows 
that any nearest neighbor VB state is a ground state. One can also mention the 2D Shastry- 
Sutherland Heisenberg model where a particular nearest neighbor VB is the unique 

ground state, and which has an experimental realization in SrCu2(B03)2 jsj. 

So far, we do not yet have any gapped liquid state.^ To obtain a qualitative idea of how 
VB states can be the building blocks of a gapped QSL, we will briefly explain the short range 
resonating valence bond (RVB) picture proposed by Anderson |g2|. If we exclude the toy 
models discussed above, a VB state is generally not an eigenstate of the Heisenberg model. 
Starting from a nearest neighbor VB state, the Heisenberg Hamiltonian will induce some 
dynamics among the VB states. If we take the kagome example, a nearest neighbor VB 
state inevitably contains some "defect" triangles without any singletj^ While the term Hijk 
leaves the VB state unchanged if the corresponding triangle has a singlet bond, the three VB 
touching i, j and k will be moved by Hijk if (ijk) is a defect triangle. The ground state 
can be viewed as a linear combination of (many) VB configurations (not necessarily nearest 
neighbor). Anderson suggested that, with appropriate interactions and lattice geometry, the 
ground state wave function could be "delocalized" over a large part of the subspace spanned 
by short range VB states. By forming a linear superposition of a large number of very 
different VB states, the system may restore all the lattice symmetries (which are broken by 
an individual VB state) and form a QSL. 

A more formal approach to this idea will be discussed in Sec. ^, but this picture can 
already be used to anticipate the nature of the magnetic excitations in such a short range 
RVB liquid. To this end, we first consider a 2D model where one ground state is equal to (or 
dominated by) one particular VB state. Contrary to the Anderson's RVB liquid, the wave 
function is localized in the vicinity of one particular VB state. It can be thought as a 2D 
analog of the Majumdar Gosh chain, where the ground state is a spatially regular arrangement 
of singlet bonds. Many 2D models are known to realize such VB crystals (VBC) p3], and 
we refer to Ref. for a recent example where the exact ground states are known. In 

a VBC, a finite energy excitation can be created by replacing a singlet bond by a triplet 
{S = 1), with an energy cost proportional to J. But is it possible to construct two separated 
spin-^ excitations in such a system ? As a trial state, one can place two remote spins "up" 
(two spinous excitations) at sites and i. Then, to minimize the energy, the regular VB 
structure of the ground state should be reconstructed as much as possible. However, due 
to the spinous, the regular pattern cannot be fully reconstructed between and i, and a 
"string" of misaligned VB is unavoidable. The unpaired spins behave as a topological defect 
in the crystalline order. So, two remote spinous perturb the ordered VB background, not 
only in their vicinity, but all the way between them. They lead to an energy cost which is 
proportional to their separation.^ So, isolated spinous are not finite energy excitations in a 

There exists a general method for constructing an SU{2) symrnetric spin model with short range interactions 
such that all the nearest neighbor VB states are ground states |28|. Building on this idea, it was possible_to 
construct SU{2) symmetric spin-^ models (with short ranged interaction) with a gapped QSL ground state ||29|| . 
Although complicated, these models are among the very few examples where the ground state is well established 
to be a short ranged resonating VB liquid. 

^"The ground states of Hmg spontaneously break the translation symmetry. On the Husimi cactus, the ground 
state is highly degenerate. The Shastry-Sutherland ground state does not break any symmetry (the ground state 
is unique), but the lattice has an even number of spins per unit cell and should be considered as a band insulator 
in our classification. 

^^Whatever the nearest neighbor VB state, exactly 1/4 of the triangles have no singlet bonds. 

^^The situation is very difl^erent in ID. In the Majumdar-Gosh model, one can get a finite energy state with two 
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VBC. The ordered VB background is a medium which confines the spinous in pairs. Since an 
RVB state should instead be viewed as a liquid (no broken symmetry, no long range order), 
it is reasonable to expect the spinons to be able to propagate as independent particles. As 
we will see in the next section, the proper way to address this question of confinement and 
deconfinement of spinons is to understand the emergence of gauge degrees of freedom in these 
systems. 



6 Schwinger bosons, large- limit, and Z2 topological 
phase 

6.1 Schwinger bosons representation 

The spin wave approach is a large-S' approach and is unable to capture highly quantum states 
which are rotationally symmetric, such as RVB wave functions. From the discussion of Sec. ^ 
it is natural to look for a description in terms of singlet fields leaving on bonds, and able 
to describe the presence or absence of a singlet between two sites. Such variables appear 
naturally when using the Schwinger boson representation of the spin operators |^, 

At each site, two types of bosons carrying a spin "up" and "down" are introduced: aj| and 

operators 



oIt'^jT ~ . St=al-faii , 5"^ = al^aif (25) 



With these relations, the commutation relations [5*^, Sf] = ie""^^ Sf are automatically verified. 
The total spin reads = + l) , where Ui — al^aif +aj^ai| is the total number of bosons 

at site i. To fix the length of the spins, the following constraint must therefore be imposed 
on physical states: 

aj^ttii + al^aii = 2S (26) 

With this representation,^ the Heisenberg interaction is of degree four in the boson op- 
erators and can be written 

Si ■ Sj = S — -{AijY Aij (27) 
with Aij — ai^Uji — aiiUj]. (28) 

The bond operators A]^ behave as a singlet creation operators: A^, when applied onto the 
boson vacuum, creates a spin singlet | TUi) ^ I iilj) and, from Eq. ^ A\^Aij is proportional 
to the number (0 or 1) of a singlet between sites i and j. In addition, Aij is invariant under 

rotations: redefining the bosons by an SU{2) matrix P: 

unchanged]^ 





^ P 











6.2 Mean field approximation 



Arovas and Auerbach | p5| suggested an approximation in which the interaction is decoupled 
using mean-field expectation values 



and to replace the constraint (Eq. 



(29) 

by a condition on the average number of boson per site 

{al^a^-f +al^a^l} = 2S. (30) 



remote spinons by introducing a domain wall in the dimerization pattern in and i. 
^^Fermions can also be used, leading to other very interesting theories for (gapped of gapless) QSL jj^, [is], [l2[ . 

" -1 



Aij can be written using the 2x2 antisymmetric tensor e - 



1 







Co-o-'iio-ajCT'- The 



rotation invariance of Aij follows from the fact that any P £ SU{2) satisfies P^eP - 
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By this replacement, the Hamiltonian becomes quadratic in the boson operator 

-"^X^i {a\-^ai^ + a\yaii - 2S\ + est. (31) 



A chemical potential has been introduced at each site to tune the boson densities so that 
they satisfy Eq. The mean field Hamiltonian Hmf (and thus its ground state |0)) depends 
on the complex parameters Qlj (one for each pair of sites ij where Jij / 0). These parameters 
have to be adjusted to satisfy the self-consistency conditions on each bond 

Q% = ^Jij{0\aii:aji - anaj^lO). (32) 

As in the spin wave approach, the Heisenberg model has been reduced to a quadratic boson 
model (here with some self consistency conditions). However, the crucial difference is that the 
present formalism does not impose any preferred spin direction: giving a finite expectation 
value A'^j 7^ to the operator Aij does not break the SU{2) symmetry - which is a necessary 
condition to describe a QSL. 

Generally speaking, two family of solutions can be found at this mean field level. In the 
first class, favored when S is large, the Schwinger boson Bose-condense in some particular 
mode. Because they carry a spin index, such condensate state (spontaneously) breaks the 
SU (2) symmetry. These solutions describe Neel states with long range spin-spin correlations. 
In such cases, the Schwinger boson mean-field theory is essentially equivalent to the spin wave 
approach (Sec. 

The second class corresponds to (mean field) QSL states. There, the ground state is 
rotationally invariant, and the Bogoliubov quasi particles obtained by diagonalizing J^mf are 
gapped. Since the corresponding creation operators, fe|,a and &i,c«, are linear combinations 
of the original bosons, these excitations also carry a spin ^. The most important question 
is whether the existence of these deconfined (free in the mean field approximation) spinons 
is an artifact of the mean field approximation, or if they could survive in some Heisenberg 
spin model. In the first a case, the inclusion of the fiuctuations that were neglected would 
confine the spinons and would deeply change the nature of the ground state. The mean-field 
picture of a fully symmetric state with non interacting spinons excitation is then qualitatively 
incorrect. Another possibility is that the spinons remains deconfined, even in presence of 
fiuctuations. In that case, the mean-field approximation is a very useful starting point. We 



will discuss in Sec. |6j a scenario where it is the case. But before, we need to introduce the 
basic formalism that is needed to describe the fiuctuations about the mean field solution, and 
emergence of gauge degrees of freedom in the system. The central question concerning the 
long distance and low energy properties of the system will be whether these gauge degrees of 
freedom confine or not the spinons. 

6.3 Large J\f, saddle point 

To discuss the role of the fiuctuations neglected in Eq. [i^, it is necessary to formulate the 
mean field approximation as a saddle point approximation in path integral formulation of the 
model. It will then be possible to identify the structure of the most important fiuctuations 
about the saddle point. To do so, one duplicates J\f times the two species of bosons (j and [). 
In addition to the site and up/down indices a, the boson operators now carry an additional 
"fiavor" index m = 1, • • • , A/". The Hamiltonian and the constraint are then generalized to 



H 



WTr E ■^^^ (33) 



2AA 

AT 



and 



m — 1 



^1; +al.a,mi = 2J^S. (35) 
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For A/" = 1, this model is Heisenberg model with SU(2) symmetry. For A/" > 1, this model 
has an enlarged symmetry given by the group 5'p(A/') S is a parameter of the model, and 
is no longer related to a representation of SU{2) if A/" > 1. The bond operator Aij is a sum 
over all the flavors. For this reason, in the limit where Af is very large, the fluctuations of Au 
become negligible compare to its expectation value and the approximation made in Eq. Ea 
becomes exact. 

A formal way to establish this result is to adopt a formulation of model where the partition 
function Z — Tr [e~^^^ at temperature T = f3~^ is expressed as a coherent state path integral 
over complex variables Zima{T) (in correspondence with the boson operators aima) which are 
periodic functions of the imaginary time r G [0, In this formalism the partition function 



reads 



26 



Z = J V[z,ma{T),X^{T)]exp(^- J^^ Lq dr^ (36) 

i m a ij 

^ ^ Xi {_Zim^ Zim,-] + ZiralZiml ~ 25') (3"^) 



1 m 



Aij — {zim^Zjml ZimlZjm^ 

) , (38) 

m — 1 

where a Lagrange multiplier A has been introduced at each lattice site and each time step 
to enforce the constraint (Eq. |5|) exactly (to simplify the notations, the r dependence of all 
fields is implicit). 

Now, a Hubbard-Stratonovich transformation is performed : 



Z = / I?[2™<,(r),A.(T),Q,j(T)]exp - / Li dr (39) 







Ll — ^ ^ ZimfrdrZima ^ ^ f j \Q'^j\ QijAij QjjA, 

Jij 



+i ^ Ai (zim| + z^mi Zimi~2S) (40) 

i m 

This new formulation involves an additional complex field Qij on each bond. The equivalence 
of Li with the initial Lagrangian Lo can simply be checked by performing the Gaussian 
integrations over Qijij) for each bond and each time step: J 2?[(5ij(r)] exp Li dr^ = 

exp ^— Lo dr^ (up to a multiplicative constant). At this point, the Af fiavors of particles 
are no longer coupled to each other, but are coupled to a common bond field Qij. So, for 
a fixed space-time configuration of Q, we have A/" independent copies of the same boson 



^® The simplectic group of 2Af x 2Af matrices Sp{Af) is the set of matrices P which satisfies P^J'P = J', where 
1 



-1 



generalizes the antisymmetric e tensor. 



1 
-1 

For an introduction to the path integral formalism in this context of quantum magnetism, see for in- 
stance Ref. pa]. We sketch the main steps of the derivation in the case of a single bosonic mode [a, at] = 1. 
For any complex number z, a coherent state \z) = e^"*|0) is defined. These states satisfy: a\z) = z\z), 
{z\z') = e^^ and the resolution of the identity / d?z \z){z\e~^''^ = 1. On writes the partition function 
as a product over Nt imaginary time steps .E = Tr ^^-d^H ^-cItH . . .j _ Xr [(1 — dTH){l — drH) ■ ■ ■] 

with dr = P/Nr. Then, the identity is inserted at each step: Z = lim]v^^oo / ^H^i '^^•^^t^ e~l^i 1^ (zi |1 — 

driT|2]v^)e~l^"x1^(2^^|l - driT|2]v^_i)---e-l^2p(22|l _ Next, we write e"!^'!^ (z^ 1 1 - drH|zi_i> ~ 

exp [—2i(zi — 2i_i) — (iTH(zi, Zi_i)], where the complex number H(z,z') = {z'\H\z) is obtained by writing the 
Hamiltonian in a normal-ordered form an replacing by z and a by z'. Taking the continuous time limit dr — > 
is formally written as Zi — Zi^i dTz{T)dT and finally leads to .2 = J 'D[z\eK'p{— Ldr) with the Lagrangian 

L = z{T)drZiT) + H{z{t), z(t)). 
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system. In addition, the Lagrangian Li is now quadratic in the z variable. We note Gq\ 
the corresponding quadratic form, a big matrix which has space (i), time (r), spin (a) and 
complex conjugacy (z versus z) indices (but no flavor index), and depends on the auxiliary 
field Q and A. L\ is then 



+ ^[2,;,(r);2,,(r)]G 



Zja' (t') 



(41) 



Performing the Gaussian integral over the z fields is now simple, as it gives (det[G])-^, also 
equivalent to e^''^''"^''^" . The partition function is now expressed as a path integral with the 
fields Q and A only, but with a complicated non-Gaussian weight: 



Z = 



0[2;,„(r),A,(r),Q,,(r)]exp \-M j La dr 



L2 = |Q,,f -2iS^A. +Tr[log(GQ„ 



(42) 
(43) 



Here, the flavor indices m have disappeared and M only appears a global multiplicative factor 
in the action. With this formulation of the Sp{M) "spin" model, it is clear that, in the limit 
A/" — > 00 the partition function will be dominated by the configurations ((5°, A") which are 
saddle points of the action 5[Q, A] = L2 dr. In other words the fluctuations of Qij and Xi 
are frozen when A/" ^ 00. Such saddle points are obtained by requiring 



dS 



dS 



= , , =0. (44) 

(45) 

and in most cases they are found to be time independent Q°j(r),Ai(r) — > Q^j,\1. The 
equations above can then be shown to be equivalent to the self consistency conditions of 
Eqs. |3^ and|32[ with = ^ I],„{0!aimTajmi - aimiajmilO). 



6.4 Fluctuations about a saddle point and gauge invariance 

We are now ready to discuss the fiuctuations that are present when M is finite, where the 
field Qij{r) is able to fluctuate around its mean fleld value Qij. Treating all the possible 
fluctuations is certainly very difficult, as it would amount to solve the original spin problem. 
A possible approach is to compute perturbatively the first 1/J\f corrections to the mean 
field results |35). However, this can miss some important effect (instabilities) which are not 
perturbative in l/A/", and will generally not shed light on the issue of spinon confinement 
that we are interested in. Instead, as in |^, we will examine the qualitative structure of 
the fluctuation modes which are important for the long distance properties of the system. In 
particular, we would like to know if some fluctuations could confine the spinous (in which case 
the mean field picture is incorrect), or if the QSL state is stable at finite A/". As we will see, 
there are some fluctuations modes which are described by a gauge fleld ji^, ^] and mediate 
some (possibly long ranged) interaction between the spinon. The dynamics of this gauge field 
is therefore crucial to the physics of the spin system. In some cases this gauge field will be 
in a confining phase, and the A/" = cxd limit (where the fluctuations are frozen out) does not 
represent the physics of the finite A/" models js^]. In some other situations, the gauge field 
has a deconfined phase and a QSL state with elementary spinon excitation is possible po||. 

First, it should be noticed that the description of the spin operators with Schwinger bosons 
is redundant in the sense that an arbitrary local change of phase in the boson operators does 
not change the physical spin operators. In the path integral formulation, this becomes a full 
space-time gauge invariance. The Lagrangian Li (Eq. Ml) is invariant under 

2ima(r) — > e^'^^^^Zimair) (46) 

Qij{r) ^ e''^'(-'+^^(-»Q,,(r) (47) 
Ht) — ^ A,(r) - c»,A,(t) (48) 
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where Ai(r) is some arbitrary angle at each site and time step. 

However, this local U{1) gauge invariance is broken to a smaller invariance group in the 
vicinity of a saddle point (Q", A"). This can be illustrated the simpler context of a classical 
ferromagnetic Heisenberg model. A ground state is magnetized in one particular direction 
and thus breaks the 0(3) symmetry of the Hamiltonian. The theory for the (transverse) spin 
deviations around this ferromagnetic state has an 0(2) symmetry, an not 0(3). The situation 
is similar for the fluctuations of the bond fleld Qij. Although the model has a local U{1) 
gauge invariance, the action describing the fluctuations around Qij have a lower invariance 
group. In the ferromagnet example, we look at the rotations under which the ground state 
is unchanged. Similarly, we look for the gauge transformations which leave unchanged. 
These transformations form the invariant gauge group (IGG) of the saddle point, a concept 
introduced by X. G. Wen A gauge transformation i Ai belongs to the IGG of Qij if 
it is static and satisfles 

Q?, = Q°,e^<^'+^^-' (49) 

If the lattice made of the bonds where Q'ij is non zero is bipartite, it is easy to show that 
Ai = 6^ on sublattice A and Ai = —6 on sublattice B satisfies Eq. ^ for any (global) angle 6. 
In such a case, the IGG is isomorphic to (7(1). On the other hand, if the lattice of the bonds 
where (5°^- / is not bipartite, the IGG is isomorphic to Z2, since Ai = tt and Ai = are the 
only two solutions to EgL ^ when Q°j / 0. 

The general result |kl|| is that, among the fiuctuations around the saddle point Q", some 
modes are described by a gauge field, with a gauge group given by the IGG. We will illustrate 
this result in the simple case IGG= Z2.pj 



6.5 Z2 gauge field 

If the IGG is Z2, the important fluctuations turn out to be fluctuations of the sign of Qij. 
We therefore parametrize these fluctuations in the following way 

Q..(r) = 0?, e'^-M , A,(r)G{0,^}. (50) 

where the field Aij will play the role of a "discrete" (Z2) vector potential living on the links 
of the lattice (pairs of sites where Qij / 0) . 

Doing the integration over all the other fiuctuation modes (amplitude fiuctuations the 
bond field Qij, fiuctuations of Ai, etc.) in order to obtain an effective action for Aij and 
the bosons Zio- onlyp] is formally possible, but it is of course a very difficult task in practice. 
One can instead determine the symmetry constraints, and, in a Landau-Ginzburg type of 
approach, construct the simplest action compatible with these symmetries. 

For this, we consider the (static) local gauge transformation i 1— > Ai with the restriction 
Ai G {0, tt}. Because Aij is defined modulo 2n, —Aj is equivalent to +Aj and the transfor- 
mation rules take the usual form (except for the discrete nature of Aij): 

Zia — > e^'-Zia (51) 
A^j — y Aj+Ai-Aj. (52) 

These local transformations form a very large symmetry group (2 to the power of the number 
of lattice sites) and severely constrain the effective Hamiltonian for these degrees of freedom. 
Because of this invariance, a term like Aij, Aij or even cos{Aij) cannot appear as an en- 
ergy termj^ Instead, only the products of e'^^'J on closed loops are gauge invariant. As a 
circulation of the a vector potential, these loop terms are the analog of the magnetic fiux in 
electromagnetism. Such products can thus appear in an effective description of the fiuctua- 
tions about the mean field solution. Terms like £ij = drAij + Xi — Xj, which are equivalent 
to the electric field, are also gauge invariant. As for the couplings to the bosons, the coupling 
to A allowed by the gauge invariance (an spin-rotations) are of the type Zi^ e^'^'^ Zja-. 

^'^The cases where IGG= U{1) are generically unstable saddle^oints: the gauge fluctuations lead to spinon 
confinement, and lattice symmetry breaking (VBC) when S = ^ ||39|. This will not be discussed here. 
^^From now on, we go back to A/" = 1 a drop the flavor index m for simplicity. 

^^In the same way, a a mass term like the square of the vector potential A^^^ is forbidden by gauge invariance in 
conventional electromagnetism. 
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6.6 A simple efFective model 



We can combine the gauge invariant terms above into a simple Hamiltonian which can phe- 
nomenologically, when IGG=Z2, describe the gauge fluctuations about a saddle point and 
their efi^ect on the spinous: 



(53) 



The operator crfj has eigenvalues ±1, like a pseudo spin-i, and corresponds to e'"^'J in the 
path integral formulation (Eq. |5^). afj corresponds to the electric field operator. In the path 
integral, Aij and £ij are conjugated. So afj and crfj should not commute on the same bond. 
The natural choice in our discrete case is f^o-fj = —cTijafj. So, afj and afj are the x and 
z components of the pseudo spin-|. The bosons represent the Bogoliubov quasi particles 
(spinon) of the mean field Hamiltonian. The first term (K) is a sum over all the elementary 
plaquettes (square here for simplicity) and corresponds to the magnetic energy of the gauge 
field. The second term (F) is the electric energy, which generates fiuctuations in the magnetic 
flux. The third one (t) describes the spinon hopping and their interaction with the gauge 
field. The last terms represents the energy cost A > to create a spinon (related to the spin 
gap of the spin model) and some (large) penalty V when more than one spinon are on the 
same site. 

This model is of course not directly related to the original spin model but contains the same 
two important ingredients that have been identifled in the large A/" limit (spinon coupled to 
Z2 gauge fleld fluctuations) and can provide as a simplifled and phenomenological description 
to a gapped QSL. 

Because of the gauge symmetry, the physical Hilbert space of the model should be con- 
strained to avoid spurious degrees of freedom: two states which differ by a gauge transfor- 
mation correspond to a single physical state and should not appear twice in the spectrum. 
In the Hamiltonian formulation of gauge theories, the solution is to construct the operators 
Uify which generate the local gauge transformations, and impose that all the physical states 
should be invariant under these transformations: Uig\phys.) — jphys.) Viq. In the present 
case, an elementary gauge transformation at site io changes the value af^j for all neighbors 
j of io (noted j G +). In addition, it changes the sign of the boson operators in io. This 
transformation is implemented by the following unitary operator 

C/.o = exp [i7r(6];^^6,o| + feIo|fe»oi)] H '^'"oi (5^) 

je+ 

The constraint Uig = 1 is the lattice version of the Gauss law, divE = p, in electromagnetism, 
and the spinous appear to play the role of the "electric" charges. 

Readers familiar with lattice gauge theories will have recognized the Hamiltonian formula- 
tion of a Z2 gauge theory p^ . However, to show that the ground state of this model realizes 
a topological phase (when F is small enough), we will show that it is very close to the toric 
code model introduced by Kitaev |H]. 



6.7 Toric code limit 

One goal of these notes was to show that (gapped) QSL in Mott insulators are topologically 
ordered states with emerging gauge degrees of freedom. To conclude, we will now take advan- 
tage of Kitaev's lectures of topological states of matter (in this school), and show the close 
connection between the large Af description of gapped QSL and Kitaev's toric code |^. 

We consider the limit of Eq. |53| when t = 0, F = and V = 00. In this limit, the bosons 
cannot hop any more, and can only be zero or one per site: rii = bl^^bigi + bl^^bigi € {0, 1}. 
Using Ui = 1 (Eq. |5^) we flnd: e"^"» = T[je+ '^fj^ so that the boson occupation numbers are 
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expressed in terms of the (lattice divergence of the) electric field operators: 2ni = 1 — njg + ""fj- 
Taking The Hamiltonian can then be written as 

H = -K'^aija^kf^liafi - ^A'^Y^'^ij (55) 

□ i je + 

which is exactly the (solvable) toric code Hamiltonian |^. 

We can now import some results from the toric code analysis. Although simple derive in 
the framework of Eq. they are highly non trivial from the point of view of the original 
spin model. First, the ground state breaks no symmetry and the spinous (here at the sites i 
with nj(=+ ""fj = ~1) frse particles, they are not confined by the gauge field fluctuations. 
Secondly, the ground state is degenerate on a cylinder or on a torus (periodic boundary 
conditions), as required by the LSMH theorem. The ground state are topologically ordered 
in the sense that no local observable can distinguish the different ground states. Beyond 
the spinous, the model also have Z2-vortex excitations, which correspond to plaquettes with 
CijO'jk^ki^H = These gapped excitations are singlet states in the original spin model 

since the bond fleld Qij and its sign fluctuations a^j are rotationally invariantPj These 
excitation have a non trivial mutual statistics with respect to the spinon, and a bound state 
of a spinon and a vison behaves as a fermion. Finally, the topological properties of the model 
(fractional excitations, topological degeneracy) are robust to perturbations, and should persist 
in presence of a small F and small t (Eq. [ssl ). 
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